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ABSTRACT 



HR 8799 is an hFO mA5 7 Doradus, A Bootis, Vega-type star best known for hosting 
four directly imaged candidate planetary companions. Using the CHARA Array inter- 
ferometer, we measure HR 8799's limb-darkened angular diameter to be 0.342ib0.008 
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mas; this is the smallest interferometrically measured stellar diameter to date, with an 
error of only 2%. By combining our measurement with the star's parallax and photom- 
etry from the literature, we greatly improve upon previous estimates of its fundamental 
parameters, including stellar radius (1.44ib0.06 Rq), effective temperature (7193ib87 
K, consistent with FO), luminosity (5.05ib0.29 Lq), and the extent of the habitable 
zone (1.62 AU to 3.32 AU). These improved stellar properties permit much more pre- 
cise comparisons with stellar evolutionary models, from which a mass and age can be 
determined, once the metallicity of the star is known. Considering the observational 
properties of other A Bootis stars and the indirect evidence for youth of HR 8799, we 
argue that the internal abundance, and what we refer to as the effective abundance, 
is most likely near-solar. Finally, using the Yonsei-Yale evolutionary models with uni- 
formly scaled solar-like abundances, we estimate HR 8799's mass and age considering 
two possibilities: 1.516lo:o24 Mq and 33l]'3 2 Myr if the star is contracting toward the 
zero age main-sequence or 1.513^q'q24 Mq and 90^*^50^ Myr if it is expanding from it. 
This improved estimate of HR 8799's age with realistic uncertainties provides the best 
constraints to date on the masses of its orbiting companions, and strongly suggests they 
are indeed planets. They nevertheless all appear to orbit well outside the habitable zone 
of this young star. 

Subject headings: planetary systems, stars: fundamental parameters, stars: individual 
(HR 8799), techniques: high angular resolution, techniques: interferometric 



1. Introduction 



HR 8799 (HD 218396, HIP 114189) is a 7 Doradus-type star, a class characterized by pulsations 
on time-scale s longer than the typical J Scu t i stars, their neighbors on the H ertzsprung-Russell (H- 
R) diagram dBreeerlEgTgl : IZerbi et allll999l : Ixave et al.lll999l : lHandlerlll9"99l l. HR 8799 is also a A 
Bootis-type star, which is a spectroscopically-defined group of non-magnetic, chemically peculiar. 
Population I A- and early F-type stars that show metal defici encies in Fe-peak elements while show- 



ing s o lar to slightly over s olar abundances of C, N, O, and S (jVenn &: Lambertlll99Cll : lPaunzen et al 



1998 



Gray &: Kaydll999l ). In a ddition, HR 8799 is a "Vega-like" star, po ssessing excess infrared 



emission longward of ~ 20 fim (jSadakane &: Nishidal Il986l : ISu et al.l l2009l ) , attributed to thermal 
emission from a debris disk. 

Aside from these three distinguishing characteristics, interest in HR 8799 was amplified in 
2008 wi th the discovery o f three planet-like companions at projected separations of 24, 38 and 
68 AU (jMarois et al.ll2008l ). From comparisons of the observed fluxes to the predictions of evolu- 
tionary models, Marois et al. estimated companion masses of 5-11 Afjupitor for companion b and 
7-13 Mjupiter for companions c and d, each. Because the masses of the companions depend on the 
age of the system, Marois et al. used four lines of reasoning to estimate the host star's age: (1) 
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comparing the star's galactic space motion to other young stars in the solar neighborhood, (2) plac- 
ing the star on a color-magnitude diagram, (3) characterizing the typical ages of A Boo and 7 Dor 
stars, and (4) the assumption that HR 8799 has a massive debris disk, which is typically only found 
associated with young stars (< 500 Myr). Based on these characteristics, they constrained an age 
range between 30 and 160 Myr. A fourth companion was imaged two years later with an estimated 



mass of ^Jupiter if the system is 30 Myr old and 10ib3 Mjupiter if 60 Myr old (IMarois et al 



2OIOI). 



Mar lev et al.l (|2012l ) fit atmospheric and evolution models to the data to determine the masses, 
radii, temperatures, gravities, and cloud properties for the planets. They found masses of 26 Mjupiter 
for planet b and 8-11 Mjupiter for planets c and d, though they acknowledge that their model fit for 
b's mass is not likely to be the true value. They also determined an age range of 360 Myr, 40-100 
My, and 30-100 Myr for each planet, respectively. 



The ages adopted by iMarois et al.l (|2008l . I2OIOI ) have become a topic of some debate. The 
star's age is of vital interest because it is directly linked to the masses of the companions. If the 
star is young, its companions are brighter and their inferred masses are lower - i.e., planets - 
than if the star is older, in which case its more massi ve companions - i.e . , bro wn dwarfs - have 
cooled significantly and could be mistaken for planets (IMoro-Martm et al.ll2010l ). For example, if 
HR 8799 is as old as the Hyades (625 Myr. iPerryman et al.lll998l ). all companions would be more 
massive than 13.6 Mjupiter j and thus more appropriately called brown dwarfs according to popular 
convention. One exain ple of the effect of presumed ages on masses involved the announcement by 
Lafreniere et al.l (|2010l ) of a directly imaged planetary companion to IRXS J160929. 1-210524 in 
the Upper Scorpius association, which relied on an assumed cluster age of 5 Myr. Soon afterward, 
Pecaut et al.l ( 20121 ) determined an older age for the cluster: 11^2 Myr. In this case, the imaged 
companion has a mass of ~14 Mjupiter, putting it just above the planetary mass limit. 



Moya et al.l (l2010al ) refute the IMarois et al.l (j2008l ) young age of HR 8799 on all four counts. 
They point out that (1) space motions of young disk stars are often inconclusive, (2) the uncertain 
internal metallicity of HR 8799 makes comparisons with evolutionary models and other stellar 
clusters unreliable, (3) that A Boo and 7 Dor type stars have a broad range of ages, and (4) that 
debris disks are "h ighly chaotic and unpredictable", and cannot be used to estimate stellar age. 



Mova et al 



Zerbi et al 



( 2010al ) instead determined HR 8799's age by mod eling 7 Dor pulsations dete c ted b y 



19991 ) and adopting spectroscopic parameters bv lAllende Prieto &: LambertI (jl999l ). 



Although primarily limited by the unknown inclination of the star's rotation axis, they find a 
preferred age of between 1.1 and 1.5 Gyr; again, this age implies that the companions to HR 8799 
are brown dwarfs. 



Moro-Martm et al.l (j2010l ) subsequently challenged the claim made by iMova et al.l (j2010a 



finding an age of the system of ~100 Myr after considering the limiting case for orbital stability; 
at these relatively small separations, massive companions would not be dynamically stable on long 
timescales. They used the double 4:2:1 mean motion resonance configuration presented by other 
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auth ors (e.g.. lReidemeister et alJl2009l : lGozdziewski &: Migaszewskill2009l : iFabrvckv k, Murrav-Clav 
20 id ) as a way to ensure dynamical stability by avoiding close encounters. Moro-Martm et al. used 
evolutionary models to predict ages from 150 to 350 Myr and conclude that the dynamical state of 
the system favors an age of ~100 Myr, which places the companions just shy of the brown dwarf 
regi me. Howev e r, the y po int out that t he ave rage temperature of 7350 K and luminosity of 5 Lq 
from ISadakanel (j2006l ) and Moya et al. ( 2010b ) implies an age for the star that is very young when 
compared to other A Boo stars. 

In addition to these recent attempts to estimate HR 8799's age, many other efforts have been 
made to determine its age in the last decade. The resulting values vary from 20 Myr up to 1.5 
Gyr. In Tabled] we summarize both the attempts described above and other efforts. A variety of 
methods were utilized, including inspecting the star's proper motion, comparing its position on an 
H-R diagram with isochrones, inspecting the star for infrared excess, asteroseismology, astrometry, 
spectroscopy, dynamical stability analysis, disk symmetry, direct imaging, and group membership. 

Here we seek to greatly improve our understanding of HR 8799's age by using interferometric 
observations to directly measure its angular diameter. This value, when combined with the HIP- 
PARCOS parallax, yields the star's physical radius. In combination with HR 8799's bolometric 
flux determined from broad-band photometry, the stellar luminosity and effective temperature are 
precisely determined. Using these results we provide new age and mass estimates based on com- 
parisons with the stellar evolutionary models, adopting a metallicity most appropriate for A Boo 
stars and HR 8799 in particular. Specifically, Section 2 describes the interferometric observations 
and calibrator star selection; Section 3 discusses the visibility measurements and how stellar pa- 
rameters were calculated, including angular diameter, radius, luminosity, and temperature; Section 
4 explores the physical implications of the interferometric observations; and Section 5 summarizes 
our findings. 



2. Interferometric Observations 



Observations were obtained using the Center for High Angular Resolution Astronom y (CHARA) 

Array, a six element optical-infrared interferometer located on Mount Wilson, California (Iten Brummelaar et al 
20051 ). All observations used the Precision Astronomical Visible Observations (PAVO) beam com- 
biner with a spectral ban dpass of ^650-800 nm . For a description of the PAVO instrument and data 
reduction procedure, see llreland et al.l (j2008l ). We observed HR 8799 over seven nights spanning 
two years with the S2-W2 and Sl-El telescope pairs with maximum baselines of 177 m and 331 m, 
respectively0 



^The three arms of the CHARA Array are denoted by their cardinal directions: "S", "E", and "W" are south, 
east, and west, respectively. Each arm bears two telescopes, numbered "1" for the telescope farthest from the beam 
combining laboratory and "2" for the telescope closer to the lab. The "baseline" is the distance between the telescopes. 
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We interleaved calibrator star and HR 8799 observations so that every target was flanked by 
calibrator observations made as close in time as possible, which allowed us to convert instrumental 
target and calibrator visibilities to calibrated visibilities for the target. Calibrators are stars with 
predicted diameters that are significantly smaller and in close proximity in the sky to the target 
star. Reliable calibrators were chosen to be single stars with angular diameter estimates <0.3 mas 
so they were nearly unresolved on the baseline used, which meant uncertainties in each calibrator's 
diameter did not affect the target's diameter calculation as much as if the calibrator star had a 
larger angular size. 

In or der to estimate the redde ning of each calibrator star, we obtained th e Tycho (Bt — Vt) 
color from iPerrvman &: ESAl (j 19971 ) and converted to {B — V) using the table in iBesselll (|200d ). We 
then compared the obs erved (B — V) value w ith the list of intrinsic {B — V) colors as a function of 
spectral type given by Schmidt-Kaler ( 1982 § to arrive at an estimate of E{B — V), and adopted 
the reddening law described in lO'Donnelll ([199J) to de-redden the observed magnitudes. The pho- 
tometric angular diameters were t hen d etermined using the relationship between the {V — K) color 
and log ^LD from iKervella et al.l (120041). Th e error in the diameters is due to the relative calibra- 



tion error stated in the 



Kervella et al 



(|2004l ) paper as well as errors in photometry measurements. 
Table [2] lists the input photometry and resulting photometrically estimated angular diameters of 
the calibrator stars. 



3. Results 



3.1. Angular Diameter Measurement 



The observed quantity of an interferometer is fringe contrast or "visibility", more formally 
defined as the correlation of two wave-fronts whose amplitude is the visibility squared (V^), which 
is fit to a model of a uniformly-illuminated disk (UD) that represents the face of the star. Diameter 
fits to V"^ were based upon the UD approximation given by V'^ = ([2Ji(a;)]/x)^, where Ji is the first- 
order Bessel function and x = ttBOudX^^ , where B is the projected baseline at the star's position, 
^UD is the apparent UD angular d iameter of the star, and A is the effective wavelength of the 
observation (jShao &: Colavitalll992l ). A more realistic model of a star's disk involves limb-darkening 
(LD), and the relationship incorporating the linear LD coefficient fix (jHanburv Brown et al.lll974l ) 
is: 



1 



2 3 



(1-^a)^^ + /^a 



1/2 J3/2(x) 
x3/2 



Table [3] lists the date of observation, the calibrator used, X/B, the calibrated visibilities (F^) 
errors in V'^ {aV'^). 



(1) 



and 



^The table was obtained from |http: / /obswww.unige.ch / gcpd/ mk01bv.html 
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The LD coefficient was obtained from IClaret &: BloemenI (l201lh after adopting the e ffective 
temperature (Teff) and surface gravity (log g) values from lAUende Prieto LambertI (| 19991 ). which 
were 7586 K and 4.35, respectively. The resulting UD and LD angular diameters are listed in Table 
[6l Figures [1] through [8] show the LD diameter fits for HR 8799 by night and by calibrator and 
Figure [9] shows all the data combined. 

Table m lists the resulting angular diameter fits for each night using each calibrator, and Figure 
[To] shows a graphical version of the results. Though there is some scatter in the diameter fit from 
each night/calibrator combination, the scatter from the 2011 data is less pronounced than from 
2010 data, which is seen in the diameter fits shown in Figures H] through [71 In particular, for 
observations obtained in 2010, the data exhibit sinusoidal-like variations about the best angular 
diameter fit; the variations are not seen in 2011 data. This is almost certainly due to coating 
asymmetries between CHARA Array telescopes (particularly overcoated silver versus aluminum) 
that were present in 2010 but removed in 2011. Attempts were made to search for polarization 
effects that would cause these residuals, but no conclusive evidence was found. Visibilities were, 
however, never measured in full Stokes parameters in a configuration that showed these residuals. 

Because the 2011 data show none of the sinusoidal residuals, our final ^ld incorporates the 
data from 2011 only. This yielded an angular diameter of 0.341ib0.008 mas. When the data from 
all nights and using all calibrators are fit together, ^ld is 0.347ib0.007 mas, which is within the 
uncertainty of the adopted value. 



The uncertainty for 0ld was calculated as described in the supplemental material from [Derekas et al 



(j201ll ): for each one of 10,000 Monte-Carlo simulated Gaussian distributions, we calibrated the in- 
strumental visibilities and fit an angular diameter to the calibrated data using a least-squares 
minimization. We accounted for random errors by adding random numbers generated from the 
empirical covariance matrix scaled by the reduced of the original fit, and then repeated the 
procedure. The final uncertainty was the resulting standard deviation of the total distribution. 

We observed HR 8799 with multiple calibrator stars every night except one to check on the 
behavior of the calibrators themselves. For example, when three calibrator stars were used, we used 
calibrator 1 to determine the angular diameters of calibrator 2 and calibrator 3 to make sure the 
stars were reliable. The calibrator HD 214698 shows the largest scatter, which is expected because 
it is the smallest calibrator of the three. There will naturally be more uncertainty when measuring 
its calibrated visibilities using the other two larger stars. The scatter seen in the night-to-night 
measurements is expected because the stars are small and very nearly unresolved. 
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3.2. Stellar Radius, Luminosity and Effective Temperature 



At a distance of 39.4 ± 1.1 pc (|van Leeuwenl l2007l l. HR 8799's Old of 0.342 ± 0.008 mas 
corresponds to a stellar radiu^ of 1.44±0.06 Rq, corresponding to a precision of 4%. We note that 
this radius is 8% large r than the radius inferred from photometric and temperat ure considera tions 



m 



2010a 



Gray &: Kavd (jl999l ) and subsequently adopted in recent asteroseismic analyses iMoya et al.l (e.g. 



b|); only 1% o f this discrepancy can be attributed to the pre-revised Hipparcos distance used 



in the calculation of lGrav &: Kavd (|l999l ). 



In order to determine the luminosity (L) and T^q of HR 8799, we first constructed its photo- 
metric energy distribution (PED) from the averages (when multiple measurements were available) 
of Johns on UBV magnitudes, S tromgren uvh y magnitudes, a nd 2MASS JiJK magnitudes, as pub- 



lished in 



mm. 



Crawford et al 



Mendoza et al 



Cutri et all (120031 ). and [Olsenl (|198.1 



(ll966l l .lBregei] (ll968ll.lEggen 



1993 



( 199dl. lHauck k Mermilliod (119981 ). ICezari et al.1 (|l999l ^. iHog et al.l (|200d l 



(ll968l l .lMermilliodl (119861) . ISchuster fc Nissen 



19941 ). Table [5] summarizes the adopted magnitudes; 



the assigned uncertainties for the 2MASS infrared measurements are as reported, and for optical 
measurements are standard deviations of multiple measurements. The assigned uncertainties of 
the optical measurements ther efore account for the low amplitude optical variability (±0.02 mag) 
observed for this variable star (jZerbi et al.lll999l ). 



The bolometric flux (-Fbol) of HR 8799 was determined by finding the best f it (via min i- 
mization) stellar spectral template from the flux-calibrated stellar spectral atlas of iPicklesI (|l998l ). 
This best PEE) fit allows for extinction, using the wavelength-dependent reddening relations of 
Cardelli et al.l (119891 ). The best fit was found using a FO V template with an assigned temperature of 



-1 —2 

cm 



7211±90 K, an extinction of Ay = 0.00 ± 0.01 mag, and a Fbol of 1.043±0.012 x lO"'^ erg s^ 
(a 1.1% precision). Figure [TT] shows the best fit and the results are listed in Table [6l To check for 
possible systematic biases in our adopted prescription, we also fit the PED using synthetic Kurucz 
spectral model^, assuming no extinction. The resulting -Fbol estimate is consistent to within the 
1.1% error reported above, further validating our technique and measured -Frol- 

The bolometric flux was then combined with the distance to HR 8799 to estimate its luminosity 
(L = Att d^FeoL), which yielded a value of 5.05 ± 0.29 Lq. The uncertainty in the luminosity (6%) 
is predominantly set by the uncertainty in the distance. The -Fbol was also combined with the 
star's ^LD to determine its effective temperature by inverting the relation. 



Fbol 



1 



^LDi^^efr) 



(2) 



where a is the Stefan-Bolzmann constant. This produces an effective temperature of 7203 it 87 K, 
determined to a precision of 1%. Because //ld is chosen using a given Tefj, we used our new Tgff 



We define 1 solar radius to be 6.960 x 10 cm, consistent with an average of previous measurements ijEmilio et al 



20121 1 



Available to download at http: / /kurucz. cfa.harvard.edu] 
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value to select ^ld and iterated, //ld increased by only 0.01 to 0.49±0.02, 0ld increased by 0.001 
mas to 0.342zb0.008 mas, and T^g decreased by 10 K to 7193 K. The very slight change in 0ld did 
not affect the radius calculation at all. 

We note that this Teg is nearly identical to the assigned temperature of the best fit stellar 
templates used to calculate -Fbol- However, our interferometric Tes is significantly less model 
dependent, and it is important to verify the accuracy of the T^g predicted by the P EP. Moreover, 
ourmeasured T^g is also consistent with the detailed spectral type classification bv lGrav fc Kave 
(|l999l ) of kA 5 hFO mA5 v A Boo, follow ing the spectral-type temperature scale for dwarf stars 
assembled in lKraus &: Hillenbrand! (|2007l ). The hydrogen lines are better tracer of the stellar T^s 
(hFO; 7200 K); the metal lines suggest a higher T^g (mA5; 8200 K) only because its atmosphere is 
slightly metal depleted. 

A po tential bias i n the size measurement of any A-type star is oblateness caused by rapid 
rotation. I Rover et al.l (j2007l ) measured a t;sini of 49 km s~^ and assuming the planets orbit in 
the same plane as stellar rotation, the actual ro tation velocity will inc rease from 49 to 104 km s~^ 
when the 28° inclination is taken into account (jSoummer et al.ll201ll ). For an A5 V star with an 
approximate mass of 2.1 Mq (ICoxll200Clll. our measured radiu s of 1.44 Rq, and the relation between 
M, oblateness, and vsini described in van Belle et al. (j2006l ). the predicted oblateness of HR 8799 
is only ~2% and thus within the 1-a errors in 0ld- Using the FO V spectral type, which is a closer 
fit to our measured Tes, the estimated mass is 1.6 Mq and the predicted oblateness is 2.6%, still a 
small effect. 



Discussion 



4.1. Habitable Real-Estate of HR 8799 



HR 8799 is currently the only directly imaged multiple planet system and there may be smaller, 
more Earth-like planets orbiting the star or even moons orbiting the imaged planets that have not 
yet been detected. This would present the possibility of life if the planets were in the habitable 
zone (HZ) of the star, so we used our new precise measurements to improve the estimate of the 
system's HZ. We used the following equations from I Jones et all (2006): 



eS) 



(4.190 X 10"^ T^g) - (2.139 x 10"^ T^s) + 1.296 



(3) 



and 



SbATes) = (6.190 X 10"9 Tjj) - (1.319 x 10"^ T,s) + 0.2341 (4) 

where Sb^iiTcs) and Sb^oiTes) are the critical fluxes at the inner and outer boundaries in units of 
the solar constant. The inner and outer physical boundaries rj o in AU were then calculated using 



/_L/Le 
Sb,iiTc 



© 



and 



Sb,oiTes) 



(5) 
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These equations assume the HZ is dependent on distance only and do not take other effects, such 
as tidal heating, into account. The inner boundary is the limit where a runaway greenhouse 
effect would evaporate any surface water while the outer bou ndary is the limit where a cloudless 
atmosphere could maintain a surface temperature of 273 K. iJones et alj (|2006l ) note that these 
equations produce a conservatively small HZ and the actual HZ may be wider. 

We obtained habitable zone boundaries of 1.62 AU a nd 3.32 AU. HR 8799's planets have 
semimajor axes of 14.5 to 68 AU (jMarois et al.ll2008l . 120101 ). There is no chance the planets orbit 
anywhere near the habitable zone unless they are highly eccentric, which is a configuration more 
likely to be unstable. 



4.2. EfTective Abundance of the A Boo Stars, and HR 8799 

The accurately determined stellar properties of HR 8799 also allow for detailed comparisons 
with stellar evolutionary models from which a mass and age can be inferred. However, correctly 
interpreting masses and especially ages depends critically upon knowing the internal abundances of 
this peculiar abundance star. As noted in the Introduction, the atmospheric abundances of Fe-peak 
elements are distinctly subsolar ([Fe/H] ~ —0.4 dex) while the abundances of C, N, O, and S are 
essentially solar. Because these surface abundances may not directly trace the internal abundances, 
it is unclear which, if any, of the available uniformly scaled abundance models to adopt for these 
comparisons; we subsequently refer to the abundance of the uniformly scaled model that predicts 
properties most consistent with observational constraints as the effective abundance. 

As an example of the effect on the inferred ages, comparisons of HR 8799 with models that 
assume subsolar abundances throughout yield ages of 1.05ib0.26 Gyr for [Fe/H]= —0.27 or 1.71ib0.18 
Gyr for [Fe/H] = —0.43, while those that assume near-solar abundances ([Fe/H] = +0.05) yield an 
age of < 0.1 Gyr; the specific physical assumptions used in these models are described in Section 
4.3. Obviously, the polyabundic atmospheres of A Boo stars and their uncertain internal abundance 
compromises the validity of these comparisons and the inferred values. Nevertheless, we present 
the available observational evidence for constraining the metallicity of A Boo stars, and HR 8799 
in particular, and suggest the effective abundance of HR 8799, and perhaps all A Boo stars, is near 
solar. 

The observat i onally favored theory to explain the A Boo phenomenon, originally proposed by 



Venn LambertI (jl990l ). is that the atmospheres have been polluted by the accretion of Fe-peak 



depleted gas. Depletion is believed to occur because o f grain format ion; a similar depletion of Fe- 



like elements has been observed in interstellar clouds (lMortonlll974l ). The grains that themselves 
accrete Fe-elements are consequently inhibited from accreting on to the star because of the stronger 
radiation pressure they experience. It is not clear why this accretion onto the grains occurs, or 
where the accreting mater ial comes from. It may be interstellar, but in many cases appears to be 



circumstellar. As noted by lGray &: Corballyl (120021 ) . all four A Boo stars within 40 pc (including HR 
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8799) exhibit an infrared excess that is interpreted as the presence of circumstellar dust, while only 
~ 18% of non-A Boo A-type stars exhibit such excesses. The apparent depletion of these elements 
requires only rela tively low accretion rates (10~^^ Mq yr~^), because of the shallow convective zones 



in A-type stars (jCharbonneaul Il99ll ). Accretion at these rates will quickly establish abundance 
anomalies within a few Myr, b ut these anomalies will likewise disappear in as little as 1 Myr once 
the accretion has terminated (jTurcotte &: Charbonneaulll993l ): the interesting implication is that 
all A Boo stars are either currently accreting, or have only recently terminated their accretion. If 
this favored theory is correct, it suggests that the depleted Fe-peak abundances do not represent 
the internal abundances of these stars. Given this and the typical limiting main-sequence lifetime of 
~2 Gyr for these intermediate mass Population I starJl, one would expect their internal abundance 
to be on average close to solar. 

We note that an alter native mechanism to explain the A Boo phenomenon proposed by 
Michaud &: Charlandl (jl986l ) suggests that the depletion of Fe-peak elements is a result of dif- 
fusion and mass loss. In this case, A Boo stars are much closer to the end of their main-sequence 
lifetimes. To be effective, the star would need to be losing mass at a rate of order 10"^'^ Mq yr~^ 
for a few times 10^ yr and thus implies that the A Boo phenomenon be restricted to the end of the 
main-sequence evolutionary phase for A stars. This predicted timescale is difficult to reconcile with 
the discovery of very young A Boo stars in the Orion OB I associations (age < 1 Myr), and the 
lack of any A Boo stars in intermediate age open clusters (jGrav &: Corballvl 12002 ). Evidence such 
as the higher fraction of these stars with circumstellar debris disks, noted above, instead suggest 
that the A Boo phenomenon is more common among young A stars. If true, this would strengthen 
the case that the internal abundances of these Population I stars, at least on average, should be 
close to solar. In particular, we note that the existence of A Boo stars in the Orion star forming 
regi on implies that their prim ordial abundances are solar, consistent with other stars in this clus- 
ter (iNieva &: Przybillal l2012l ). If the depletion of Fe-like elements is restricted to only the outer 
atmospheres, then one can conclude that in these cases the internal abundances should likewise be 
solar. 



As described in lMarois et al.l ( 20081 ). the space motions of HR 8799 are consistent with those 
of young stars in the solar neighborhood. We investigated this further by c omparing the UVW 
space motions of HR 8799 to those of nearby moving groups assembled in iTorres et al.l (|2008l ): 
kinematic association with a moving group would not only help constrain the age of the system 
but also the primordial abundance of the star. We calcula ted UVW space motions using HR 8799's 
HIPPAR COS distance and parallax (Ivan LeeuwenI 120071 ) . and a radial velocity of —12. Oil. 3 km 
s~ 
V 



1 from iBarbier-Brossat FigonI (|200d ). This yielded velocities of f7 = -12.24±0.37 km s-\ 
= -21.22±1.10 km s~^ and W = -7.15±0.86 km s"^; these are all within '-I km s"^ of the 



values reported in lMarois et al.l (|2008l ). Figure [12] illustrates the space motion of HR 8799 relative 
to that of several kinematically similar moving groups after adopting moving group velocities and 



'''Studies have shown that stars currently forming in open clusters are generally metal-rich (e.g., 



Santos et allboogi ). 
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uncertainties from [Torres et al.l (|2008l ). Although the space motion of HR 8799 is not consistent 
to within la of any of these groups, it is consistent to with in 3a of two of thes e groups: Columba 
{1.2a), and e Cha {2.2a). We thus confirm the report by iMarois et alj (|2008l ) of the potentially 
youthful kinematics of HR 8799, which favors but can not confirm an age less than < 1 Gyr, and 
likely a solar abundance consistent with many of these stars. 

As described in the Introduction, there are two other aspects of HR 8799 that argue for a young 
age for the star and system. The str ong infrared excess of HR 8799 statistically favors a n age less 
than ~ 500 Myr (jGaspar et al.ll2009l ). despite the challenge raised bv lMova et al.l (12010^1 . An even 



strong er youthful age restriction is imposed by the dynamical considerations of iMoro-Martm et al 



(|2010l ). Although their conclusions require use of the predictions of poorly constrained planetary 
evolutionary models, even considerable uncertainties in the mass estimates require the age of the 
system to be less than a few 100 Myr, otherwise the planets would be too massive to be orbitally 
stable. The ages much less than 1 Gyr for HR 8799 can only be reconciled if the adopted effective 
abundances are near solar. 



Recently, iMova et al.l (|2010bl ) used the 7 Dor-type pulsations of HR 8799 to asteroseismically 
constrain the star's internal metallicity. They find a best fit internal metallicity of —0.32 it 0.1, 
assuming that the rotation axis of the star, which affects the rotation speed of the star, is somewhat 
highly inclined relative to our line of sight (50°). This hi gh inclination is inconsistent with the 
more face-on orientation of HR 879 9's disk (ISu et al.ll2009l ) and the apparent orbital plane of its 
companions (jLafreniere et al.ll2009l ). although these need not be coplanar. The asteroseismology 
analysis is compromised by the use of an assumed radius that is too small by 8% and a temperature 
that is too hot by 3% relative to the precisely determined values in this study. The measured R and 
Tgff presented here will be valuable to constrain the model parameter space in future asteroseismic 
modeling efforts, which will possibly constrain the internal metallicity. 



Moya et al.l (|2010bl ) do find possible - though less likely - solutions with internal abundances 
close to solar ([Fe/H] = —0.12) for closer to pole-on orientations; this is the most metal-rich 
metallicity reported in their study. Overall we find the results of this metallicity analysis to be 
inconclusive primarily because of the large uncertainty in the inclination of the star's rotation axis, 
which Moya et al. acknowledge is the limiting factor in their analysis. If anything, the result 
that they find some acceptable solutions with near solar metallicities and young ages corroborated 
previous asteroseimo logy studies of A Boo stars. Using densities inferred from stellar pulsations, 
Paunzen et al.l (|l998l ) found the location of most A Boo stars showing 6 Scuti pulsations on a plot 
of average density versus period to be consistent with the positions of 5 Scuti stars with solar 
abundances. They interpret this as evidence that the low Fe-peak abundances are restricted to the 
surface of A Boo stars and do not represent the state of the interior composition. Although the 
variations in HR 8799 are co mpatible with 7 Dor pulsations, the large number of discovered hybrid 
6 Scu - 7 Dor pulsators (e.g.. lRowe et al.ll2006l : iGrigahcene et al.ll2010l ) suggests that this argument 
could be applicable for all A Boo stars. 
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In summary, based on the observational evidence supporting the accretion of clean gas theory 
explaining the abundance pattern of A Boo stars, the existence of at least some A Boo stars in an OB 
Association, the evidence of youth for many A Boo stars (and especially HR 8799), the restriction of 
ages less than ~2 Gyr for Population I stars, and the available asteroseismic constraints on internal 
metallicity, we conclude that the most appropriate effective abundances for HR 8799, and quite 
possibly all A Boo stars, is near solar. Until more sophisticated evolutionary models that account 
for the polyabundic atmospheres and possibly interiors of these stars are developed, we are hopeful 
that A Boo stars in clusters or with lower mass companions will be discovered, which will enable 
improvements in both assigning effective abundances and tests of the validity of using uniformly 
scaled abundances. 



4.3. Mass and Age of HR 8799 and Implications for Its Companions 



Following the above arguments, we estimated t he mass and ag e of HR 8799 using the stellar 



evolut ionary models of the Yonsei-Yale group (Y^, lYi et al. 



2001 



updated by iDemarque et al 



(|200J) to account for convective core overshoot. The models handle convection using mixing length 
prescription, adopting a mixing length of 1.7432 times the pressure scale height that is set by 
comparisons with a solar mod el. The model uses the solar abundances of iGrevesse et al.l (119961 ) 
and OPAL radiative opacit i es (Rogers Sz Iglesias 1995 : Iglesias &: Rogers 19961 ) for the interior and 
the [Alexander &: Fergusonl (|1994| ) opacities for the cooler, outer regions of the s tar. Add i tiona l 
info rmation regarding the in put physics assumed in these models can be found in lYi et al.l (|200ll ) 
and IDemarque et al.l (|2004l ). However, since a set solar metallicity models is not provided, we 
followed the recommendation of the modelers themselves and generated a set of solar metallicity 
models by linearly interpolating isochrones and mass tracks between nearest the metallicity models 
of [Fe/H]=-0.27 dex and [Fe/H] = +0.05 dex (models x74z01 and x71z04 models, respectively). 
The solar metallicity models are illustrated in Figure [131 

At HR 8799's temperature, stars are predicted to have fully contracted to their smallest size, or 
the zero age main sequence (ZAMS) at 40 Myr, and then begin expanding due to stellar evolutionary 
effects at a much slower rate afterwards. As such, HR 8799 is either contracting onto the ZAMS 
or expanding from it, and we use the 40 Myr isochrone as the upper limit for the age in the pre- 
main-sequence scenario and as the lower limit for the age in the post-main-sequence scenario. More 
specifically, we find that if HR 8799 is contracting onto the ZAMS, it has an age of 331(3 2 Myr 
and a mass of 1.516^g Q24 -^o- If HR 8799 is expanding from the ZAMS, we find it to have an age 



of 90tfo^ Myr and a mass of 1.513 



+0.023 
-0.024 



Mq. These masses and ages are found by interpolating 
betwee n the solar rnetallic ity models described above. The masses are consistent with the mass 
used in lMarois et al.l (|2008l ). which was 1.5+0.3 Mq. We remind the reader that these quoted errors 
on the star's mass and age are statistical and therefore do not take into account uncertainties in 
the metallicity and/or models themselves. 



With an age of < 0.1 Gyr, the companions that HR 8799 harbors are more likely to have 
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planetary masses. As explained in iMarois et al.l (|2008l ). planetary mass objects could only be as 
bright as these observed companions are if they are young. That inferred age, and thus companion 
masses, depends critically upon adopted evolutionary model abundance. While near solar seems 
likely, even slightly subsolar abundances can give uncertainties in age that extend to above ~ 500 
Myr, increasing the inferred companion masses by at least a factor of two. 



5. Summary 

We measured the angular diameter of HR 8799 using the CHARA Array interferometer and 
used our new value of 0.342ib0.008 mas to calculate the star's physical radius (1.44ib0.06 Rq), 
luminosity (5.05ib0.29 Lq), and effective temperature (7193ibl0 K) by combining our measurement 
with information from the literature. We used our T^g measurement to determine the size of the 
habitable zone, which is well inside the orbits of any of the companions detected to date. 

Based on a variety of techniques, we concluded that the most appropriate abundances for HR 
8799 are close to solar. We combined our R and Teg values with isochrones to estimate the 
star's mass and age in two scenarios: 1.5161qq24 Mq and 33l]'3 2 if the star is contracting onto 
the ZAMS or 1.513l[]:[j^^ Mq and 90^50^ if it is expanding from it. In either case, this young age 
implies the imaged companions are planets and not brown dwarfs. The only case in which the 
companions would be close to brown dwarf mass is if the highest age of 471 Myr was the true case, 
and even then the masses would be on the exoplanet/brown dwarf cusp. 

We would like to thank Gerard van Belle for his insight on the nature of HR 8799's pirouette 
through space. The CHARA Array is funded by the National Science Foundation through NSF 
grant AST-0606958 and by Georgia State University through the College of Arts and Sciences, 
and STR acknowledges partial support by NASA grant NNH09AK731. This research has made 
use of the SIMBAD database, operated at CDS, Strasbourg, France. This publication makes use 
of data products from the Two Micron All Sky Survey, which is a joint project of the University 
of Massachusetts and the Infrared Processing and Analysis Center /California Institute of Tech- 
nology, funded by the National Aeronautics and Space Administration and the National Science 
Foundation. 
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Table 1. HR 8799 System Parameters from the Literature. 





Age 






(Me) 


(Myr) 


Method Used 


Reference 


1.47±0.30 




spectral synthesis & spectrophotometry 


Grav & Kave (1999) 




50-1128 


H-R, diagram placement 


Song et al. (2001) 




30 


stellar kinematics, H-R diagram 


Zuckerman & Sone (2004) 






placement vs. isochrones 






20-150 


Local Association membership 


Moor et al. ("2006") 




590 


IR excess 


Chen et al. (2006) 




30 


IR excess, H-R diagram placement 


Rhee et al. C2007~) 


1.5 


30-160 


multiple methods 


Marois et al. C2008') 


1.2-1.6 


<100 Myr 


dynamical stability analysis 


Gozdzicwski & Mieaszewski (^2009") 




^50 


observational data analysis 


Rcidcmeister et al. (2009) 


1.5 


<100 Myr 


dynamical stability analysis 


Fabrvckv & Murrav-Clav ('2010) 




~100 Myr 


disk inclination 


Moro-Martfn et al. (2010") 


1.32-1.33'' 


1123-1623 


asteroseismology 


Mova et al. f2010a) 


1.44-1.45'' 


26-430 


asteroseismology 


Mova et al. (■2010a) 


1.32 


1126-1486 


asteroseismology 


Mova et al. (2010a") 




30 


Columba Association membership 


Dovon et al. (2010") 




30 


Columba Association membership 


Zuckerman et al. (2011) 




30-100 


direct imaging 


Currie et al. (2011) 




30-300 


atmospheric/evolution model fitting 


Marley et al. (2012) 




30-155 


dynamical stability analysis 


Sudol & HaehighiDour (2012) 



Note. — ''Solution 1; ''Solution 2 
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Table 2. Calibrator Star Properties. 





Spectral 


V 


K 


E(5 - V) 


^photometric 


HD 


Type 


(mag) 


(mag) 




(mas) 


213617 


Fl V 


6.42 


5.58±0.02 


0.03 


0.288±0.006 


214698 


A2 V 


6.33 


6.21±0.02 


0.05 


0.189±0.005 


219487 


F5 V 


6.60 


5.54±0.02 


0.05 


0.304±0.006 



Note, 
are from 



Mermil 



Ferryman &: ESA 



Spectral typ es are from SIMBAD; V magnitudes 
iodi (1 199 in except for HD 214698, which is from 
(|l997f K No errors were listed so we assigned errors 



of ±0.01 mag; K magnitudes are from lCutri et al.l (|2003l ). 
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Table 3. Calibrated Visibilities. 





Calibrator 


Spatial Preq 




Date 


Used 


(10^ rad-i) 





2010 Aug 25 HD 213617 


2, 


.223 


0, 


.861 


0. 


,081 




2 


.245 


0, 


.768 


0, 


,069 




2, 


.267 


0, 


,785 


0, 


,069 




2, 


.289 


0, 


,807 


0, 


,076 




2, 


.311 


0. 


,804 


0. 


,089 




2, 


.332 


0. 


,729 


0. 


,094 




2, 


.353 


0, 


,811 


0, 


,138 




2, 


.375 


0, 


,844 


0, 


,146 




2, 


.397 


0, 


,924 


0, 


,142 




2, 


.419 


0, 


,976 


0, 


,135 




2, 


.442 


0. 


,892 


0. 


,092 




2, 


.464 


0. 


,834 


0. 


,066 




2 


.486 


0, 


,788 


0, 


,059 




2, 


.507 


0, 


,764 


0, 


,057 




2, 


.528 


0, 


,804 


0, 


,063 




2, 


.550 


0, 


,742 


0, 


,059 




2, 


.571 


0. 


,728 


0. 


,067 




2, 


.593 


0. 


,713 


0. 


,084 




2, 


.613 


0, 


,771 


0, 


,104 




2, 


.634 


0, 


,840 


0, 


,120 




2, 


.655 


0. 


,784 


0, 


,103 




2, 


.677 


0. 


,827 


0. 


,096 




2, 


.699 


0. 


,801 


0. 


,077 



Note. — All data were taken using the S2-W2 baseline 
(177 m). Only a portion of this table is shown here to demon- 
strate its form and content. A machine-readable version of 
the full table is available on the online version of The Astro- 
physical Journal. 
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Table 4. HR 8799 Angular Diameter Measurements. 





Calibrator 


^UD 




Date 


HD 


(mas) 


(mas) 


zUiU Aug zo 


01 Qd'i T 


n on'T-Ln noi 
u.zy ( ±U.Uzl 


n Qnnj-n noi 
U.oUy±U.UZi 




OI CiA QT 


n 001 _Ln m 
U.ool±U.Ulo 


n onT-Ln m a 
O.oy7±0.01b 




A 11 

All 


U.oo/±U.UiO 


n Q'TO-Ln m 
U.o/z±U.Ulb 


zUlU oep U< 


OI Qi^l 7 


U.o4z±U.UiD 


U.O0b±U.Ui4 




OI O/l 


U.O/5±U.UlU 


n om _i_n m n 

o.oyi±o.oio 




A 11 

All 


o.3oy±u.uoy 


n o^o_i_n m n 
0.O/O±0.010 


zUlU bep Uo 


01 O^il T 

zldblT 


U.30z±U.U18 


n OI c j^n m n 
0.31o±0.019 




01 CiA OT 


n oi o_Ln m o 


n oo^j_i_n m o 
O.ozb±0.01o 




A 11 

All 


U.oUoztU.Uil 


n 0O1 _i_n m o 


om 1 on 
zUii UCt zU 


OI Qi^l 7 

ziobi / 


n Qon_i_n m q 


U.ooo±U.Ui4 




91 /1RQ8 
Z14Dyo 


n '^/iQ-i-n ni 7 
u.o4iyitu.ui I 


n QfiQ-un ni s 

U.oDoitU.UiO 






n Q 1 /) _i_n m 
U.314±U.Ul0 


n oo^_i_n m k 
O.oz r ±0.015 




A 11 
All 


n OQ1 _i_n nnn 

u.ooi±u.uuy 


f\ o A cc_i_n m n 
O.o4o±0.010 


om 1 /~\«4- oi 
zUll Uct zl 


OI O^il T 

zldblT 


n oo£j j^n ni n 
U.33b±U.UlU 


n oc;n_Ln mn 
0.350±0.010 




zi4byo 


n Qon_i_n m o 
U.ozU±U.Uiz 


n OQO_i_n m o 
0.OOO±0.01o 




91 Q/1S7 


n QQ9-i-n ni a 

U.OOZztU.UlD 


n Q/iR-i-n ni 7 




A 11 
Ail 


n Qon_i_n m n 

u.oou±u.uiy 


n O/io^n noi 
O.o4o±0.0zi 


2011 Oct 22 


213617 


0.308±0.010 


0.321±0.011 




214698 


0.308±0.014 


0.320±0.015 




219487 


0.332±0.008 


0.346±0.009 




All 


0.328±0.023 


0.341±0.022 


2011 Sep 30 


214698 


0.323±0.009 


0.338±0.008 


All dates 


213617 


0.327±0.008 


0.341±0.008 


Aug 2011 data only 




0.326±0.008 


0.340±0.008 


All dates 


214698 


0.327±0.009 


0.341±0.010 


All dates 


219487 


0.343±0.008 


0.358±0.008 


Aug 2011 data only 




0.328±0.009 


0.342±0.009 


All dates 


All 


0.333±0.007 


0.347±0.007 


Final fit: 2011 data 


All 


0.327±0.007 


0.342±0.008 
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Note. — Figure [10] shows a graphical view of these values. 
Data were obtained using the S2-W2 baseline at 177 m for all 
nights except for 2011 Sep 30, which were obtained using the Sl- 
El baseline at 331 m. 
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Table 5. HR 8799 Photometry 





Average 


No. of 


Band 


(mag) 


measurements 


U 


6.191±0.016 


2 


B 


6.235±0.016 


2 


V 


5.980±0.016 


2 


u 


7.471±0.025 


7 


V 


6.468±0.025 


7 


b 


6.143±0.027 


7 


y 


5.962±0.026 


7 


J 


5.383±0.027 


1 


H 


5.280±0.018 


1 


K 


5.240±0.018 


1 
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Table 6. HR 8799 Stellar Parameters. 



Parameter Value Reference 

V magnitude 
K magnitude 
TT (mas) 
Distance (pc) 

/fA 

0.00 ± 0.01 This work; PED fit 
i^BOL (10""^ erg s"^ cm~2) 1.043 ± 0.012 This work; PED fit 



^cff, estimated 


(K) 7211 ± 90 


This work; PED fit 


6'uD (mas) 


0.327 ± 0.008 


This work 


6'ld (mas) 


0.342 ± 0.008 


This work 


-^linear (Rq) 


1.44 ± 0.06 


This work 


Teff (K) 


7193 ± 87 


This work 


L (Lq) 


5.05 ± 0.29 


This work 




If the star is contracting towards the ZAMS: 


Age (Myr) 


00 + 7 

— 13 2 


This work 


Mass (Mq) 


1 trif5+0:038 
-L-O-LD_0.024 


This work 




If the star is expanding from the ZAMS: 


Age (Myr) 


qn+381 


This work 


Mass (Mq) 


1.513l°;°23 Mo 


This work 



5.98±0.01 
5.24 ± 0.02 
25.38 ± 0.70 
39.40 ± 1.09 

0.49±0.02 



Mermilliod (1991 1 
Cutri et al. (2003) 
van Leeuwen (2007 ) 
van Leeuwen (2007) 
Claret k Bloemen (2011) 
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Fig. 1. — ^LD fit for HR 8799 on 2010 Aug 25. The diamonds and vertical lines are the measured 
visibilities and their associated errors, the solid line is the best-fit LD diameter, the dashed line is 
the 1-a error, and the dotted line is the 2-a error. 
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2,2 



2.3 



2.4 



2.5 



2.6 



2.7 




2.1 



2.2 



2.3 



2.4 



2.7 




SPATIAL FREQUENCY (10*' rod"') 
Fig. 2. — ^LD fit for HR 8799 on 2010 Sep 07. The symbols are the same as in Figure [H 



-27- 



1 A 
1 .? 
1 .0 
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0.2 
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Calibrator HD 213617^ 



<> <> <> 
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2.7 



> 1.4 
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0.0 



:alibrator HD 21948: 



1.2 



1.3 



2.4 



2.5 



2.6 



2.7 



1 .4 
1 .2 
1 .0 
0.8 
0.6 
0.4 

0.2 
0.0 



All calibrators 



<> 



<> 
< > 



<> <> 



If <> 



<>...(>. 



2.2 



2.3 2.4 2.5 2.6 



2.7 



SPATIAL FREQUENCY (10*' rod"') 
Fig. 3. — ^LD fit for HR 8799 on 2010 Sep 08. The symbols are the same as in Figure [H 
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2.2 



2.3 



2.4 



2.5 



2.6 



> 

h- 
< 

m 

< 



0.9 

0.6 

0.3 
0.0 



2.1 



2.2 



Calibrator HD 214698: 



2.3 



2.4 



2.2 



2.4 



2.6 



2.7 




2.5 



2.6 




SPATIAL FREQUENCY (10*' rod"') 
Fig. 4. — 6'ld fit for HR 8799 on 2011 Oct 20. The symbols are the same as in Figured! 
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:alibrator HD 219487^ 



2.1 



2.2 



2.3 



2.4 



2.5 



2.6 




SPATIAL FREQUENCY (10*' rod"') 
Fig. 5. — ^LD fit for HR 8799 on 2011 Oct 21. The symbols are the same as in Figured! 
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Calibrator HD 214698: 
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Calibrator HD 219487^ 
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2.4 
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2.6 
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2.1 



2.2 
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SPATIAL FRECUENCY (10^ rod"') 
Fig. 6. — ^LD fit for HR 8799 on 2011 Oct 22. The symbols are the same as in Figured! 
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201 1 data 
All data 

All data using cols HD 213617, HD 214698 
Ca HD 213617, all data 



Co HD 213617 

Ca HD 214698 

Ca HD 219487 

Co HD 219187 



2011 data 
all data 
all data 
2011 data 



All data 



By calibrator, 
all nights 



07/25/10 
07/25/10 
09/07/10 
09/07/10 
09/08/10 
09/08/10 



ca HD 

ca HD 

CO HD 

ca HD 

ca HD 

CO HD 



21 361 7 
219487 
21 361 7 
21 9487 
21 361 7 
21 9187 



By night/ccI combo, 
2010 data 



08/20/1 1 
08/20/1 1 
08/20/1 1 
08/21 /I 1 
08/21 /I 1 
08/21 /I 1 
08/22/1 1 
08/22/1 1 
08/22/1 1 
09/30/1 1 



CO HD 

ca HD 

ca HD 

ca HD 

CO HD 

ca HD 

CO HD 

CO HD 

ca HD 

CO HD 



21 361 7 
21 4698 
219487 
21 361 7 
214698 
219487 
21 361 7 
214698 
21 9187 
214698 



0.25 



0.30 



0.35 

LD DIAMETER (mas) 



0.40 



By night/cal combo, 
2011 data 



0.45 



Fig. 10. — A comparison of ^ld fits by individual night and calibrator. The vertical solid line 
represents the final adopted 9hD (0.342 it 0.008 mas, the top line) and the vertical dashed lines are 
the errors in that fit. Table U] lists the numerical values. Note the reduced scatter in the 2011 data 
versus the 2010 data. 
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E 
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CO 
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\— 

CD 




cn 

CD 

o 



0.1 7 
0.0 - 



P -0.1 - 
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Fig. 11. — HR 879 9 PED fit. Upper panel: The solid-line spectrum is a FO V spectral template 



from IPicklea ()1998l ). The crosses indicate photometry values from the literature. The horizontal 
bars represent bandwidths of the filters used. The X-shaped symbols show the flux value of the 
spectral template integrated over the filter transmission. Lower panel: The crosses are the residuals 
around the fit in fractional flux units of photometric uncertainty. For more details, see section 4. 
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V VS. U W vs. U W vs. V 




-16 -14 -12 -10 -8 -16 -14 -12 -10 -8 -24 -22 -20 -18 -16 -14 

U (km/s) U (km/s) V (km/s) 



Fig. 12. — The UVW space velocities of HR 8799 (black point) with 1-a (error bars) and 2-cr (gray 
elhpse) errors. Also plotted are the UVW space velocities and 1-a errors (aselhpses) of four young 
steUar associations with similar kinematics, as taken from [Torres et al.l (|2008l ). HR 8799 matches 
the space velocity of the roughly 30 Myr old Columba association to 1.2c7. 
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0.25 i [Fe/H]= 0.00 



of 0.20 



° 0.15 



0.10 -I 




1.7 M© 



1.6 M 







0.02 Gyr 

0.04 Gyr 

0.20 Gyr 

1.00 Gyr 



1.4 M> 



3.90 



3.88 3.86 
log(T) in K 



3.84 



Fig. 13. — Radius and tempe rature of HR 8799 plott ed along with isochrones (blue lines) and mass 
tracks (solid red lines) from iDemarque et al.l (120041 ) . with solar abundances. Note that the mass 
tracks predict that these stars are still gravitationally settling at 0.02 Gyr, reach their smallest size 
at 0.04 Gyr, and expand in size thereafter. 



